Thick brane solutions and topology change transition on black hole backgrounds 
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We consider static, axisymmetric, thick brane solutions on higher dimensional, spherically sym- 
metric black hole backgrounds. It was found recently [lj|, that in cases when the thick brane has 
more than 2 spacelike dimensions, perturbative approaches break down around the corresponding 
thin solutions for Minkowski type topologies. This behavior is a consequence of the fact that thin 
solutions are not smooth at the axis, and for a general discussion of possible phase transitions in the 
system, one needs to use a non-perturbative approach. In the present paper we provide an exact, 
numerical solution of the problem both for black hole- and Minkowski type topologies with arbitrary 
number of brane and bulk dimensions. We also illustrate a topology change transition in the system 
for a 5-dimensional brane embedded in a 6-dimensional bulk. 
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I. INTRODUCTION 

Curvature corrections to Dirac-Nambu-Goto (DNG) 
membranes m higher dimensional black hole back- 
grounds attracts a lot of attention recently in several dif- 
ferent areas of modern physics. The simplest corrections 
to the DNG brane action originate from small thickness 
perturbations that are quadratic at the leading order in 
the thickness of the brane [f|. These type of corrections 
have been studied lately in the case of a static, axisym- 
metric brane on the background of a spherically symmet- 
ric black hole in arbitrary number of dimensions [l], || . 

This, so called, brane - black hole (BBH) system,in the 
infinitely thin case, was first introduced by Frolov [7[ as a 
toy model for the study of merger- and topology changing 
transitions between higher dimensional black type solu- 
tions (or phases) of the Einstein field equations [sl-fTol] . 
The model also turned out to be very useful through the 
AdS/CFT correspondence for the study of phase transi- 
tions in certain strongly coupled gauge theories [ill, EH , 
while the higher dimensional generalizations of the Bern- 
stein conjecture [H, [ljjj, the study of the stability of 
brane - black hole systems [lf| , and the question of possi- 
ble micro black hole formations in high energy collisions 
(like of the LHC) [l~8[ are also important directions which 
are dealing with similar setup, and provide motivation for 
the study of the curvature corrected problem. 

Stiffness (or thickness)corrections to the BBH system 
has been studied first in [6| , with a perturbative approach 
near the critical solution of the thin system in the Rindlcr 
zone. It was found that when the brane has more than 
2 spatial dimensions, supercritical solutions behave quite 
differently from subcritical ones, and no evidence for the 
existence of such solutions has been found. A possible ex- 
planation of this behavior was that stiffness corrections to 
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the brane action break the symmetry between the super- 
and subcritical solutions, and quantum-gravitational ef- 
fects might cure the problem. 

The thickness corrected BBH system has been fur- 
ther studied in [l[ within a more general perturbative 
approach that also looked away from the Rindler zone. 
In accordance with Q , it was found that no regular sub- 
critical thick solutions exist if the brane has more than 
2 spacclikc dimensions. This result however obtained a 
simple explanation with observing the fact that pertur- 
bative approaches break down around the corresponding 
thin solutions at the axis of the system, because the zero 
order (thin) solutions are not smooth there. As a conse- 
quence, the topology change transition in the thick sys- 
tem can not be studied in the general case within a per- 
turbative approach, and one needs to find a new, exact 
solution of the problem for Minkowski type topologies. 

In the present paper we provide an exact numerical so- 
lution of the thick-BBH problem with arbitrary number 
of brane and bulk dimensions for both Minkowski- and 
black hole topologies. This work is an organic continua- 
tion of our previous perturbative approach and we kindly 
refer the reader to [if for those definitions, notation and 
results that might be missing here and would make the 
present paper completely self-contained. 

The plan of the paper is as follows. In Sec. II we 
run trough the same quick overview of the infinitely thin 
case that we provided in [l|, and reintroduce the most 
important parts of the BBH setup to make the paper 
self-contained. In Sec. Ill we write up the curvature cor- 
rected DNG action and introduce the correction param- 
eters. In Sec. IV we derive the Eulcr-Lagrangc equation 
and the equation of motion for the thick BBH system. 
In Sec. V the regularity conditions are discussed while in 
Sec. VI the far distance asymptotic solution is obtained. 
In Sec. VII we provide the exact, numerical solution of 
the problem in the near horizon region for both topolo- 
gies, and in Sec. VIII a topology change transition is 
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analyzed by considering the energy properties of a quasi- 
static thick brane evolution in the case of a 5-dimensional 
brane in a 6-dimensional bulk. 



where, and throughout this paper, a dot denotes the 
derivative with respect to r, and the action (0]) reduces 
to 



II. THE THIN BBH SYSTEM SETUP 

Let us overview, in this section, the important proper- 
ties of the BBH system, introduced in Q , that we intend 
to study in the presence of a small brane thickness in 
the following sections. We consider static brane config- 
urations in the background of a static, spherically sym- 
metric bulk black hole. The metric of an TV-dimensional, 
spherically symmetric black hole spacetime is 



ds 2 = g ab dx a dx b = -fdt 2 + f^dr 2 



r 2 dn%_ 2 , (1) 

where / = f(r) and dil%_ 2 is the metric of an TV — 
2 dimensional unit sphere. One can define coordinates 
8i{i = 1, . . . , TV — 2) on this sphere with the relation 



<m 2 i+1 = do 2 



7 <+ 



idn 2 



(2) 



The explicit form of / is not important, it is only assumed 
that / is zero at the horizon ro, and it grows monotoni- 
cally to 1 at the spatial infinity r — > oo, where it has the 
asymptotic form [19| . 



1 = 1- 



ro 



N-3 



(3) 



In the zero thickness case, the test brane configura- 
tions, in an external gravitational field, can be obtained 
by solving the equation of motion coming from the Dirac- 
Nambu-Goto action 



S = J d^-defry^ , 
where 7 is the induced metric on the brane 



lav = 9ab- 



dx a dx b 



(4) 



(5) 



and (^(n = 0, . . . , D — 1) are coordinates on the brane 
world sheet. The brane tension does not enter into the 
brane equations, thus for simplicity it can be put equal 
to 1. It is also assumed that the brane is static and 
spherically symmetric, and its surface is chosen to obey 
the equations 



9n = --- = 9 



N-2 



= it/2 



(6) 



With the above symmetry properties the brane world 
sheet can be defined by the function On-i = 0{r) and we 
shall use coordinates C M on the brane as 



e = {t,r,(j> 1 ,...,(j} n } with n = D-2 



(7) 



With this parametrization the induced metric on the 
brane is 



l„ v dCdC = -fdt 2 + 



1 2 

7 + r 



dr 



r 2 sin 2 (M7 2 , 



(8) 



S = AtA n / Co dr , 



C = r n sin n 6\J 1 + fr 2 6 , 



(9) 



(10) 



where At is the interval of time and A n = 2n n / 2 /T(n/2) 
is the surface area of a unit n-dimcnsional sphere. 



III. THE THICK BRANE ACTION 

In the case of a thick brane, the curvature corrected 
effective brane action is obtained by Carter and Gregory 
Q, and the corrections to the thin DNG action are in- 
duced by small thickness perturbations 



S= j d D ^-det 1 ^[-%{l + C 1 R + C 2 K 2 )\, (11) 

where R is the Ricci scalar, K is the extrinsic curvature 
scalar of the brane and the coefficients C\ and C2 arc 
expressed by the wall thickness parameter £ as 



Cl -^'^ 2 



c» = -li- 



te 3 

The parameter (i is related to the thickness by 



t 



1 



^V2A 



(12) 



(13) 



which originates from a field theoretical domain-wall 
model, where /i is the mass parameter and A is the cou- 
pling constant of the scalar field. 

After integrating out the spherical symmetric part and 
the time dependence on the introduced static, spherically 
symmetric, higher dimensional black hole background, 
one obtains (see also [l|) 



S = AtA n / C dr 



C = -^Co[l + s8] , 
where we introduced the notations 



£ = T^, S = aK 2 + bQ, 



(14) 
(15) 

(16) 



with 



Q = K a h K h a , a = — —L 2 , b = - — —L 2 , (17) 

b a; 24 24 V ^ 



and L is the relevant dynamical length scale of the system 
which has to be large compared to the thickness param- 
eter I. The explicit expressions of the curvature scalars 
K and Q are given in (35) and (36) of [l|. 
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IV. THE EULER-LAGRANGE EQUATION 

The curvature corrected DNG-brane action has a de- 
pendence on the second derivative of 9, thus the Eulcr- 
Lagrangc equation of the problem reads 



d 2 


( dC ) 


d 


fdC\ 


8C 


dr 2 


\de) 


dr 


Kdb) 


' 89 



(18) 







Plugging the thickness corrected Lagrangian density ([15 
into (fTS"]) we get 







d_ fd£o\ _ dCo 
dr\ 89 J 89 

^_ fa(c sy 

dr 2 \ 



d 
dr 



d(C S)\ , d(C S) 



89 J dr \ 89 J 09 
from which the actual equation of motion becomes 



(19) 



3(4) 



+ Ti0< 3 > + T 2 (6, 9, 9, , f\ /, /, r) = , (20) 



where 



1 



rfF 2 



(2n + 4)/ + Arf + 2nrf cot ( 



+ r 2 f 2(n- 
+ 2w 3 / 2 cot 



3)/ - rf 



- 10r i f 2 9 } 



(21) 



and T2 is given in the appendix. 

The goal of this paper is to provide a regular, ex- 
act, numerical solution of (|20j) in arbitrary brane and 
bulk dimensions for both Minkowski- and black hole type 
topologies, and study the properties of a quasi-static 
brane transition between the two type of solutions. 



V. REGULARITY CONDITIONS 

A. Minkowski embedding case 

As a first step we consider the asymptotic behavior 
of (|20|) near the axis of the system (9 = 0) in the 
Minkowski embedding case, where the perturbative ap- 
proaches failed to provide regular solutions in the dimen- 
sions n > 2. The asymptotic equation (obtained by tak- 
ing the series expansion of (f2T)|) around 9 = 0) reads 



£3 

t>3 



£2 
/)2 



Si 



so 



, 



(22) 



where the functions si, S2, and S3 are given in the ap- 
pendix. For a regular solution one needs to require S3, 
S2 and si to disappear at r\. 

The requirement S3 | ri = is automatically satisfied in 
the cases of n = and n = 2. For other dimensions we 
obtain the formula 



— (6 + an) 



|ri ± V 2(2a + b-an)r 2 f 



(23) 



We are looking for real solutions, hence we have the con- 
dition 



2(2a + b - an) 



> 



(24) 



Plugging the explicit values of a and b into (|24|) one can 
find that the inequality is satisfied for all n > 3. In 
the case of n — 1, the regularity condition (f2"3")l has no 
real solution for 9\ ri , which is a rather unexpected result. 
Note that this result is independent of the explicit values 
of a and b. 

From the requirements S2\ ri = and Si| n = we 
obtain unique conditions for 9\ ri = g{9\ ri ) and 0^| ri = 
h(9\ ri , 9\ ri ), where the functions g and h are given in the 
appendix. 

Accordingly, one can always find a regular solution for 
the equation ((20)) at the axes (9\ ri — 0) if n > 3. In this 
case the initial conditions are also uniquely fixed by the 
regularity conditions. In the cases of n = and n = 2, 
regular solutions do exist, however one has a freedom in 
fixing the initial condition 0\ ri . Having fixed 9\ ril the 
remaining conditions are uniquely determined. 

In the n = 1 case, the regularity requirements can not 
be satisfied with the presented method. It is interest- 
ing however, that in this case the perturbative approach 
worked well, and provided unique, regular solution to 
the problem. Furthermore, as explicitly constructed field 
theoretical domain- wall models [ly, [l7| clearly show the 
existence of regular solutions in this case, we believe that 
the obtained irregularity has no any physical meaning, 
rather it originates from some kind of break down of the 
applied method in the special case of n = 1. A similar 
problem occurred for example during the integration of 
the far distance solution in [l[, where in the n = 3 case, 
the solution developed a resonant term. Unfortunately, 
since the equation of motion in the present exact case is 
so complicated and so highly non-linear, it is very diffi- 
cult to follow up analytically the source of this irregular 
behavior, and thus its origin is presently not clear to the 
author. 



B. Black hole embedding case 

In the black hole embedding case the perturbative ap- 
proach provided regular thick solutions in every brane 
dimensions. Nevertheless, for a complete analysis, we 
examine the asymptotic behavior of (|20|) near the black 
hole horizon and compute the exact, numerical solution 
of the problem. 

As r — > 7"o, (i.e. on the horizon), the metric function 
f(r) goes to zero (see dS}) but all of its derivatives are 
nonzero and finite. We can thus consider the asymptotic 
behavior of (|20[) by taking its scries expansion around 
f(r) = 0, and get 



P 



f 



yo 







(25) 
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where y\ and yi are given in the appendix. Similarly to 
the Minkowski embedding case, from y 2 one can obtain 
a unique condition for 9\ ra = G(6\ ro ), and another for 
6» (3) | ro = H(e\ ro ,e\ ro ) from y x . The functions G and H 
are given in the appendix. In addition one can fix 9\ ro 
arbitrarily between and 7r/2 (which we will actually do 
later to consider a quasi-static brane evolution), however 
one still remains free to fix the initial condition 9\ ro . 

Accordingly, the black hole embedding problem can 
be regularly solved in any dimension, but the regularity 
requirements do not fix the initial conditions uniquely. 
For being able to provide an exact, numerical solution, 
we will adopt the initial condition for 9\ ro from the per- 
turbative approach where it was uniquely fixed by 
regularity considerations. This choice is supported by 
the fact that in the black hole embedding case there is 
no problem with the perturbative approach around the 
thin solution, and the exact thick solution is dominated 
by the linear term in the e-expansion on the horizon. The 
contributions of the higher order terms to the 9\ ro initial 
condition are negligible. 



and so on, and look for the solution of the transformed 
asymptotic equation 



„W + 2(4^(3) _ k v „ 







(28) 



as t — > 0. In (|28p a prime denotes a derivative with 
respect to t. 

Even though we observe that (|2"8")l has an irregular sin- 
gular point at t = 0, one can easily find a completely 
regular, exact solution 



v{t) = Pt 



y{r) = - 



(29) 



This is exactly the far distance behavior that we ex- 
pected, and the coefficient P, just as in the thin case, 
can be referred as the distance of the brane from the 
equatorial plane at infinity. 



VII. NUMERICAL SOLUTION IN THE NEAR 
HORIZON REGION 



VI. THE FAR DISTANCE SOLUTION 

As r — > oo the far distance solution can be obtained 
from the condition that the brane behaves asymptotically 
as a D — 1-dimensional plane. We look for the solution 
in the form 



77 i \ 
2 



(26) 



where v(f) is small and we keep only its linear terms in 
(|2"U|) . The asymptotic equation reads 



y ffl + 2(n+2 );/ ( 3 ) _ k ..^_ fc(n+2)^ 



, (27) 



where k — (8 (a + 6)e) _1 . 

It is instructive to compare (f2"T|) to the asymptotic form 
of the perturbation equation (57) of Q. As one would 
expect, the terms up to the second derivative are iden- 
tical, while the 3rd- and 4i/i-order terms in (|2~7} arc the 
explicit correspondents of the source term that appears 
in (57) of 0. 

For the solution, we obviously expect the same be- 
havior that we obtained in the perturbative approach, 
although (|2"7| can not be integrated in a simple closed 
form. In order to check our expectation, we map ana- 
lytically the point at infinity into the origin using the 



inversion transformation method (see eg. [21|) 



After obtaining the regularity conditions and providing 
the far distance solution of the exact curvature corrected 
problem, we consider the numerical solution of (|20[) in 
the near horizon region. Since we found that the initial 
conditions are not completely fixed for Minkowski em- 
bedding topologies in the cases of n = and 2, and also 
concluded that one can not satisfy the regularity condi- 
tion (|2"4")l in the n = 1 (sheet) case with the presented 
method, we choose to illustrate the numerical solution 
in the case of a D = 5-dimcnsional (n = 3) thick brane 
embedded in a 6-dimensional bulk. 

Within the perturbative description (see [l| for details) 
we could not approach the curvature singularity arbi- 
trary close unless we adjusted the thickness of the brane 
accordingly by changing the value of the perturbation pa- 
rameter e in every step. Instead, the method we followed 
was fixing e to its maximum value, which is equivalent to 
choose the thickest possible brane configuration for a pre- 
viously obtained length scale parameter L that had been 
determined by the boundary position (9q or n) of the 
brane. As we approached the singularity, the absolute 
maximum of the perturbations were growing and finally 
we arrived to a minimal 9q value, at which the applied 
perturbation method reached its limitation. 

With the present, exact solution however, nothing pre- 
vents us to approach the singularity as close as we please, 
and the only restriction we have to bear in mind is the va- 
lidity of the curvature corrected DNG-brane action pT|) , 
i.e. 



1 

d ^ 2 d 
dr dt 

cP_ . A d 2 

dr 2 



dt 2 



dt 



< 1 



By keeping the concept of obtaining the thickest possible 
brane configuration for given boundary data, just as in 
the perturbative approach, we choose to fix the pertur- 
bation parameter e to its maximum value obtained from 
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the condition 



0.1 



L 2 



0.01 



This requirement provides us the thickest possible brane 
configurations for every given length scale determined by 
the boundary data 8 with no restriction on how close 
we are to the curvature singularity. 



A. Minkowski embedding solutions 

In the perturbative approach, regular thick brane solu- 
tions did not exist for Minkowski topologies if the brane 
had more than two spacelike dimensions, i.e. for n > 2. 
In Fig. 1, we have plotted a set of exact thick brane 
solutions with varying boundary condition (6 for black 
hole embeddings and r\ for Minkowski embeddings) in 
the case of a 5-dimensional (n = 3) brane embedded in a 
6-dimcnsional bulk. 




FIG. 1. The picture shows a sequence of subcritical and su- 
percritical thick brane solutions in the case when a D = 5 
dimensional brane embedded into a N = 6 dimensional bulk. 
The different configurations belong to different initial condi- 
tions of r\ and #o, respectively. For simplicity, the bulk black 
hole's horizon radius is put to be 1, and R and Z are the 
standard cylindrical coordinates. 

There is no apparent difference on Fig. 1 if we com- 
pare it to the corresponding thin solutions. This is sim- 
ply because we are considering second order curvature 
corrections to the thin brane action so the effects are in- 
deed very tiny. To make these effects visible we plot a 
sequence of the difference function 

A6(r)=6(r)-6 DNG (r), 

of the corresponding thick and thin exact solutions, 
which is the analog of the perturbation function tp(r) 
defined in (41) of [J. 

The qualitative behavior of the individual AO curves 
are very similar to the <p perturbations that we obtained 
in the n = 1 case in [l| . The curves that are very close to 
the horizon has a sharp negative maximum and change 
sign before decaying. On the actual shape of the curves 
however we can observe the effect of the higher orders 
as having an extra oscillatory pattern. In Fig. 2 we have 



plotted a AO curve corresponding to r% — 1.001 boundary 
condition. 



AO 




1.010 



1.015 



1.020 



FIG. 2. The picture shows a near horizon A8 curve with 
minimum horizon distance ri = 1.001. 



If we go away from the horizon, the difference function 
becomes more spread and doesn't change sign before de- 
caying. The negative maximum and the oscillatory extra 
pattern are still present. In Fig. 3 we have plotted a A9 
curve corresponding to r\ = 2 boundary condition. 



AO 



-o.oi 



-0.02 



-0.03 - 



-0.04 - 



-0.05 




2.030 



FIG. 3. The picture shows a Ad curve a bit further away 
from the black hole horizon with minimum horizon distance 
n = 2. 

As we can see, there is no essential difference in the 
individual A9 curves compared to the perturbation func- 
tion tp obtained in [l[ in the n = 1 case. A remarkable 
difference appears however, if we take a look on the abso- 
lute maximum of the curves as we increase the minimum 
distance parameter r\. In Fig. 4 we have plotted a se- 
quence of AO curves. It can be seen that the absolute 
maximum of the difference function curves does not de- 
cay monotonically with increasing distance from the hori- 
zon. Instead, it has a growing tendency until it reaches 
a maximum around r\ ~ 1.15, after which it starts de- 
caying monotonically as one would expect. The concrete 
numerical values of the absolute maximum of the curves 
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AO 



AO 




FIG. 4. The picture shows a sequence of A8 curves with 
increasing minimum distance parameter in the 1.001 < ri < 2 
interval. 



are not relevant as they depend on the chosen normal- 
ization condition, i.e. on the maximum length parameter 
L. In analyzing the Minkowski embedding solutions, we 
have used the normalization parameter that belongs to 
the boundary condition 0o = jg- 

Apart, of course, from the very existence of the thick 
solutions in the n > 2 cases, the above non-monotonic 
decaying is an interesting and clear difference between 
the exact and perturbative solutions. 



B. Black hole embedding solutions 

The non-monotonicity in the maximum values of the 
difference function A8 with respect to the initial value 
parameter 6> 6 (0, is also present in the black hole 
embedding case. As we mentioned earlier, the regularity 
requirements do not fix the 9\ ro conditions for the black 
hole case, and to provide a set of exact solutions we took 
the missing conditions from our previously obtained per- 
turbative results 

We also discussed that with the presented exact solu- 
tion we can approach the curvature singularity in prin- 
ciple as close as we please, and it is only the capacity of 
our computing facility that can put some limitation on 
us. Taking into account our limits, we have chosen to 
solve our equation starting from the parameter 9$ = J^ , 
which is a significant increase compared to the limit yg, 
what we had within the perturbative approach. The cor- 
responding AO curve is plotted on Fig. 5. 

The oscillatory pattern of the higher order effects are 
clearly present just as in the Minkowski embedding case, 
and before decaying eventually, the curve change its sign 
3 times. This can not be seen on Fig. 5 as the amplitude 
is decaying very rapidly and the second and third sign 
change takes place further away from the horizon. 



As we increase the 9q parameter from 



!)()(> 



to 



(i.e. 




-1.x 10 



FIG. 5. The A6 function in the black hole embedding case 
with 8o — at the very near horizon region. 



and the maximum of the corresponding A9 curves are 
changing. The tendency in these features is not so sim- 
ple as in the Minkowski embedding case, for example 
there are several local extrema of the maximum with re- 
spect to 6*o- Nevertheless, to illustrate some of the solu- 
tions, we picked out three examples with initial param- 
eters 9 ~ |, ^ and |. The corresponding curves are 
plotted on Figs. 6, 7 and 8, respectively. 

A8 

0.000015 r 




1.15 



1.20 



1.25 



1.30 



changing the initial position of the brane on the hori- 
zon from the near pole region to the equator) the shapes 



FIG. 6. The A6 function in the black hole embedding case 
with O ~ Jl. 

For all curves in the black hole embedding case, the 
normalization has been calculated according to the length 
scale parameter corresponds to 9 Q = 



VIII. QUASI-STATIC PHASE TRANSITION 

Having in hand the exact numerical solution of the 
curvature corrected BBH problem, we are in the position 
to analyze the topology changing transition between the 
black hole and Minkowski embedding phases by consider- 
ing the energy properties of a quasi-static brane evolution 
from an initially equatorial configuration. This method 
has been introduced first by Flachi et al. in [2(|,and its 
details for our purposes have been explained in [IJ. 

The idea is very simple, one compares the energy of a 
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-0.00001 



S 1.10 1.15 1.20 1.25 1.30 



FIG. 7. The AO function in the black hole embedding case 
with do =± . 



AO 




10 1.15 1.20 1.25 1.30 



-0.00001 



FIG. 8. The A6 function in the black hole embedding case 
with #o — . 



thick brane configuration that belongs to a chosen initial 
parameter 8o (in the black hole embedding case) or r\ 
(in the Minkowski embedding case) to the energy of the 
equatorial configuration, i.e. to the brane with 8o = \ ■ If 
one plots the energy difference AE introduced this way as 
a function of the cylindrical distance parameter defined 
as 

Zoo = cos6'(r 00 ) 

(where r M denotes a certain cut off distance very far 
from the horizon) for every initial parameter 0q and n, 
then the obtained plot exhibits a loop, or instability zone, 
that is a typical sign of a first order phase transition in 
dynamical systems. For further details please refer to [l[ 
or 0. 

In Fig. 9 we have plotted the AE(Zoa) curves for 
the presented 5-dimensional brane solutions in a 6- 
dimensional bulk. The red curve represents the evolution 
of the brane in the black hole embedding phase, while the 
blue belongs to the evolution in the Minkowski embed- 
ding phase. It can be seen on the picture that something 
interesting happens around Z m ~ 0.8, where the two 
curves seems to touch each other. 

If we enlarge this part on Fig. 10, it becomes appar- 
ent that indeed the two curves overlap in the 0.76 < 




FIG. 9. The picture shows the AE(Zoo) curves of a quasi- 
static thick brane evolution in the case of a 5-dimensional 
brane in a 6-dimensional bulk. The red curve belongs to the 
black hole embedding-, while the blue curve to the Minkowski 
embedding evolution. 



Zoo 0.77 interval, which clearly shows that there is 
an energy degeneracy in the brane evolution at the near 
singularity region. In fact this overlap exhibits a real 
one-dimensional loop in the sense that both curves has 
a turn-back point in this region. The black hole em- 
bedding (red) curve increases from zero until it reaches 
its maximum around Zoo — 0.77, where it turns back 
and start decreasing and converges to some point in the 
overlap region. The Minkowski embedding (blue) curve 
starts decreasing from the same point (this point repre- 
sents the curvature singularity) until Z M ~ 0.76, where 
it turns back and starts increasing to infinity. Unfortu- 
nately, since being a one-dimensional loop, it can not be 
seen on the graph. 



AE 

320000 



290000 




1 1 2. 

0.762 0.775 
FIG. 10. The enlargement of the overlap region on Fig. 9. 



In any case, the above property of the brane evolution 
demonstrates that the corresponding topology changing 
transition remains a first order one in the curvature cor- 
rected BBH system, just as it was in the infinitely thin 
case. 
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IX. CONCLUSIONS 

In the present paper we further studied the effects of 
higher order, curvature corrections to the dynamics of 
higher dimensional brane black hole systems. Since ear- 
lier results [l|, @ clearly showed that perturbative ap- 
proaches fail to provide regular solutions near the axis of 
the system in Minkowski type topologies, we considered 
a different, exact, numerical approach to the problem. 
We analyzed the asymptotic properties of the complete 
4i/i-order equation of motion of the thick BBH system, 
provided its asymptotic solution for far distances, and 
obtained regularity conditions in the near horizon region 
for both Minkowski and black hole embeddings. 

We showed that the requirement of regularity for 
the thick solution defines almost completely the bound- 
ary conditions for the Eulcr-Lagrange equation in the 
Minkowski embedding case. The only exceptions are the 
brane configurations with I, 2 and 3 spacelike dimen- 
sions. In the case of I and 3 spacclike dimensions regular 
solutions of the problem exist, however, which is an un- 
expected result, with 2 spaclikc dimensions the problem 
can not be solved with the applied method. Neverthe- 
less it is obvious that the lack of a regular solution in 
this case should not have any physical reason as a regu- 
lar perturbative solution clearly exist in accordance with 
explicitly constructed field theoretical domain wall solu- 
tions [la, Q3 to the problem. 

For the black hole embedding solutions, the boundary 
conditions are not completely determined form the regu- 
larity requirements, nevertheless, since the perturbative 
approach, presented in [l], is valid for this case, the miss- 



ing conditions can be borrowed from those results. 

Having in hand all the boundary conditions we pro- 
vided a set of exact, numerical solutions in the near hori- 
zon region for both Minkowski and black hole topologies 
in the case when a 5-dimensional brane embedded into 
a 6-dimcnsional bulk. Based on the obtained solutions 
we discussed the properties of a topology changing tran- 
sition in the system, and concluded that the transition 
remains first order just as it was in the thin model. 

The main results of the paper are the construction of 
the exact numerical solution for the thick BBH problem 
in any dimensions (apart from the n = 1 special case), 
and the clarification of the question whether the simplest, 
higher order, curvature corrections can modify the order 
of the phase transition in the system. 

The obtained results are relevant in several research 
areas that are considering similar setup to the BBH sys- 
tem. A few of them were mentioned in the introduction 
section. 
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Appendix: Coefficient Functions 



64(a + b)ef 2 r 4 F 4 

- 20 2 /V(24O(a + 6)0 2 0e(-3 + n)r 2 - 480(a + b)¥ W 3 + 800 2 £r(-71(a + 6) 

+ 45an + 426n + 2 (a - 26)n 2 + 24(a + 6)0nr 2 cot 6) + 40 3 e(6(-8O + (69 - 7n)n) 
+ o(-80 + n(54 + (29 - 3n)n)) + 20(a(37 - 17n) + 6(25 - 8n))nr 2 cot 0) 
+ 46 6 nr 3 cot 0(2en(6 - an) - r 2 + er(-3/(a + 2b + 3an) - 4a/r) 

- 2e{b - a(-2 + ra))(-2 + n) esc 2 ) + 40 r 2 (en(-a(-25 + n)n + 6( — 3 + lira)) 

+ (5 + 4n)r 2 + er(/(26(8 + n) + a(13 + n(27 + 5n))) + r(/(5a - 36 - 2an) - 2a/ (3) r)) 
+ en(6(13 - lira) + a(28 + (-27 + n)n)) esc 2 ) + 40 4 r(era(6(57 - 31ra) 
+ a(89 - n(26 + 19n))) cot + 0r 2 (2era(5an + 6(-l + 4n)) + r 2 + er(3/(a + 26 + 3an) 
+ 4a/r) + 2e(a(6 - 5n) + 6(5 - 4ra))racsc 2 )) + /0 7 r 5 ((3a + 6)/e(l + n)r + 2{2aen 
+ r 2 + {a + b)fer 2 - 2ae(-l + ra)ra esc 2 ))) + 2/ 2 r(3(2a + b)f 3 6 ? er 9 

- 23(a + b)f 2 d 6 enr 7 cot + 80er(a(2 + ra)(4 + 5n) + 6(8 + n(15 + 4n)) 

+ 6(a + 6)0nr 2 cot0) + 80e(n 2 (36 + a(2 + ra)) - 90(a + 6)0 2 /r 5 + 0ra(6(25 + 8n) 

+ a(28 + 17n))r 2 cot 0) - 2/0 5 r 5 (2(a + 46)en 2 + 9r 2 + 3er(/(56(-l + n) + a(4 + lira)) 

+ 8(o + b)fr) - 2(a + 46)e(-l + n)ncsc 2 ) + 40 3 r 2 (en(6(3 - 9ra) + ara(-21 + 19ra)) 

+ (-7 - 4n)r 2 + er(2r(/(6(9 + 4ra) + 2a(5 + 7ra)) + (5a + 26)/ (3 V) + /(2(6(-78 

+ ra(13 + 2ra)) + 2a(-42 + ra(-5 + 6ra))) + (47a + 266)0nr 2 cot 0)) - e(-l + ra)ra(-96 

+ a(-24 + 19ra)) esc 2 ) + 20 4 r 3 (3(7a + lO6)0/ 2 er 4 + 2erar(/(6(47 + 8n) + a(74 + 53ra)) 

+ 4(5a + 26) fr) cot + 6ra cot 0(era(-36 + an) + 2r 2 + e(36 - a(-4 + n))(-2 + ra) esc 2 )) 

+ 40 2 r(en(216n + a(-6 + n(50 + ra))) cot + 0r 2 (4(a - 26)en 2 - 3r 2 + 6er(/(a(-48 + ra) 

+ 6(-46 + 3n)) + (a + b)fr) - 4(a - 26)e(-l + n)n esc 2 ))) - 40 6 /'V 7 (20e(a + 6 

+ on) (-24 + n(17 + ra)) + 20e(-13 + 7n)(a + 6 + an)r - 20 2 £n(-13 + 7n)(a + 6 

+ an)r cot + 3 r 2 (e(l + ra)r(/(3(a + 6) + 5an) + 2a/r) + 2(2era 2 (a + 6 + an) 

+ (1 + n)r 2 - 2e(-l + n)n(a + 6 + an) esc 2 ))) - 2/(-9(2a + 36)/ 3 5 £r 8 

- 8e(6 - a(-2 + ra))(-l + ra)ra cot + 2(-17a + 46)/ 2 4 enr 6 cot - 40r 2 (-r 2 

+ 2er (,/(20a(2 + n) + 6(40 + 17n)) + (10a + 96) /V) + en(6 - 7an) cot 2 + 8aera esc 2 ) 

- 40r(en 2 (-6 + 7an) + (-2 - n)r 2 + 2er(r(2/(3a(2 + n) + 6(5 + ra)) + (2a + 6)/ (3 V) 
+ /(2a(2 + n)(2 + 3n) + 6(8 + n(15 + ra)) + (28a + 196)0nr 2 cot 0)) + e(b 

+ a(8 - 7n))n 2 esc 2 ) + 2/0 3 r 4 (en(-21an + 6(3 + 4n)) + 7r 2 + er(/(2a(62 + n) 

+ 6(102 + 7n)) + 126/r) + en(6 - 46n + 3a(-8 + 7n))csc0 2 ) - 26>V(-3(46a + 436)0/ 2 er 4 

+ 2n cot 0(-en(56 + an) + 4r 2 + er(/(74a + 416 + 4(4a + 6)n) + (5a + 46) fr) 

+ e(-2 + n)(56 + a(4 + ra)) esc 2 ))) + /V(-480(a + 6)0 2 0e(4 + n)r 2 

- 160(a + 6)0 3 er 3 + (a + 6)/ 3 9 er 9 - 2(5a + 6)/ 2 8 £ nr 7 cot 

- 800 2 er(6(24O + (47 - 16ra)n) + a(240 + (68 - 13ra)ra) + 36(a + 6)0nr 2 cot 0) 

+ 80 3 e(6(-48 + 13(-4 + n)n) + a (2 + ra)(-24 + n(-18 + 5n)) + 30(a + 6)0 2 /r 5 

+ 20n(166n + a(-9 + 34n))r 2 cot 0) + 80 4 r(era(6(-42 + 43n) + a(-78 + n(76 + 15ra))) cot 

- 0r 2 (16(a + 6)en + 3r 2 + er(f(-Ua - 236 + 43an + 166ra) + 6(3a + 2b) fr) 

+ e(6(19 - 16ra) + a(16 - 13?i))ncot0 2 )) + 20V((5a + b)8f 2 er A + 4n cot 0(-76en + 5aen 2 
+ 4r 2 + er(/(26(7 + 2ra) + 7a(l + 6n)) + (17a + 4b) fr) + e(76 + a(12 - 5ra))(-2 + n) csc0 2 )) 
+ 2/0 7 r 5 (er(/(-56(2 + n) - a(14 + 25n)) - 2(9a + 5b) fr) + 2(en(6 - llara - 46n) - 5r 2 
+ en(6(-5 + 4n) + a(-12 + lira)) csc0 2 )) - 40 5 r 2 (er(-2/(-8a + 6 + 19an + 46n)r 

- 4(4a + 6)/ (3) r 2 + /(a(76 + (150 - 19ra)ra) + 6(58 + (9 - 4ra)ra) + 6(3a + 46)0rar 2 cot 0)) 
+ 2(en(a(44 - 15ra)ra + 6(-6 + 17n)) + (9 + 6ra)r 2 + en(6(20 - 17n) 

+ a(50 + ra(-62 + 15n))) esc 2 ))) + 4(-2(10a + 76)/ 3 3 £ r 6 + / 2 er 3 (2(18a + 176)0r 

+ 0(2Oa(2 + n) + 66(6 + ra) + (26a + 76)0rar cot 0)) + 2n cot 6{-en(b + an) + r 2 - 2a/er 2 

+ e(-2 + n)(6 + an) csc0 2 ) + 2/r(2en(2a + 6 - 2an) cot0 

+ 6r(-r 2 + 4(2a + 6)/er 2 + en((6 - 7an) cot 2 + 8a esc 2 )))) ] 
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S3 = 



S2 



(a + 6)r 4 / 2 
1 

16(a + 6)r 3 / 2 



n(n - 2)(1 + r 2 /0 2 )(6 + an + 2(6 - a{n - 2))r 2 /0 2 ) 



n(2/0(6 + a{8 - 7n))r + 40 f(n - I) (a + b + an)r 4 



+ f(29n(8a + 6 - 7an) + 20(8a + b~ 7an)r + f6 {7an + 4bn - 8a - 56)r 3 ) 
+ 20 2 ./'V(20(6(4?i - 5) + a(5n - 6))r + 0(6(7n - 9) + a((27 - 5n)n - 24) 
+ a/0 2 (n-l)r 3 ))) 
1 



si 



32(a + 6)r 4 e/ 2 F 2 



n(-4en(6 + an) + 4r 2 + 40 e/ 4 (-13 + 7n)(a + b + an)r 



+ 2sr(-4afr + /(8a + 46 - 8an + (26a + 76)/0 2 r 3 )) + 40 2 / 3 r 4 (40 2 e(6(-ll + 
+ a(-19 + n(8 + 3n))) - 48(a + b)0' ''er 2 + 29er(9(-37a - 256 + 17 an + 8bn) 
+ 8(a + 6)0 (3) r) + 4 r 2 (2e?i(-6 + an) + r 2 + er(3/(a + 26 + San) + 4a fr))) 
+ f 2 r 2 (4,9 2 e(b(2 + 19n) + a(-2 + n(44 + 3n))) + 48(a + 6)0 2 £r 2 
- 12(3a + 4b)6f6 3 er 4 - (5a + 6)/ 2 6 er 6 + 89er(9(28a + 256 + 17cm + 86n) 
+ 8(a + 6)0 (3) r) + 46»V(en(-56 + 3an) + 3r 2 + er(/(46(2 + n) + a (4 + 33n)) 
+ (13a + 46)/r))) + 2/(4e(6 - a(-2 + n))(-l + n) + 0r 2 (4(28a + 196)0/er 2 
+ 9{-8bsn + 6r 2 + er(2(7a + 46) /r + /(786 + 86n + 20a(7 + 2n) 
-(9a+116)/0 2 r 3 ))))))] , 
g(e\ ri ) = -((0(2(6 + a(8 - 7n))(/n + fr) + 2$ f 2 (-l + n)r 4 (2/(a + b + an) + afr) 
+ 9 2 fr 2 (-2/(6(9 - 7n) + a(24 + n(-27 + 5n))) + /(6(-5 + 4n) 

+ a(-8 + 7n))r)))/(2fr(b + a(8 - 7n) + 20 2 /(6(-5 + 4n) 
+ a(-6 + 5n))r 2 ))) | n , 
1 



D^ki) 



[ (4en(b + an) - 4r 

64(a + b)9ef 2 r i F 2 

- 40 6 e/ 4 (-13 + 7n)(a + 6 + an)r 6 + 2er(Aafr + /(-46 + 8a(-l + n) 

- (26a + 76)/0 2 r 3 )) - 40 2 / 3 r 4 (40 2 £ (6(-ll + 6n) + a(-19 + n(8 + 3n))) 
+ 200e(6(-25 + 8n) + a(-37 + 17n))r - 48(a + 6)0 2 £ r 2 

+ 9 A r 2 (2en{-b + an) + r 2 + er(3/(a + 26 + 3cm) + 4a fr))) 

+ 2/(-4e(6 - a(-2 + n))(-l + ra) - 4(28a + 196)0/0er 4 

+ (9a + 116)/ 2 4 er 6 - 20 2 r 2 (-46e?i + 3r 2 + er(/(70a + 396 + 4(5a + b)n) 

+ (7a + 4b) fr))) + / 2 r 2 (-40 2 £ (6(2 + 19n) + o(-2 + n(44 + 3n))) 

- 800e(6(25 + 8n) + a(28 + 17n))r - 48(a + 6)0 2 £r 2 + 12(3a + 46)0/0 3 £r 4 
+ (5a + b)f 2 9 6 er 6 - 40V(en(-56 + 3an) + 3r 2 + er(/(46(2 + n) 

+ o(4 + 33n)) + (13o + 46)/r))))] | n , 



1 



2/2 



2(10a + 76)/ 3 er b + / 2 er 3 (2(18a + 176)0r + 0(2Oa(2 + n) 



16(a + b)er 4 

+ 66(6 + n) + (26a + 7b)9nr cot 0)) + 2n cot 9{r 2 - 2a/er 2 

+ e(6 + an)(-n + (-2 + n) esc 2 0)) + 2/r(2e(6 - 2a(-l + n))n cot 

+ 0r(-r 2 + 4(2a + b)fer 2 + en((6 - 7 an) cot 2 + 8a esc 2 9)))}\ r , 
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Vi = ™, \n a [ (98a + 836) / a 6»"er 8 - 2(35a + l%b) f 6* enr b cot0 (A. 
32(a + 6)er 4 

+ 46V(en 2 (-6 + 7an) + (-2 - n)r 2 + 2er(r(2/(3a(2 + n) + 6(5 + n)) 

+ (2a + 6)/ (3) r) + /(2a(2 + n)(2 + 3n) + 6(8 + n(15 + n)) 

+ (28a + \%)6nr 2 cot (9)) + e(b + a(8 - 7n))n 2 esc 2 0) 

+ 2f8 3 r 4 (en(b - 7an - 46n) - 3r 2 + er(/(-6a(34 + 7n) - 6(174 + 19n)) 

- 4(8a + 76)/r) + e?i(6(-5 + in) + a(-8 + 7n)) esc 2 0) 

+ 26> 2 r 2 (-(210a + 1976)0/ 2 er 4 + 2n cot 0(en(-36 + an) + 2r 2 

+ er(/(66a + 376 + 4(6a + b)n) + (9a + 46) /V) 

+ e(36 — a(— 4 + n))(— 2 + n) esc 6»)) + 4(32(a + 6)0 (3) / 

+ 2e(6 - a(-2 + n))(-l + n)ncot + 0r 2 (-r 2 + 2er(/(20a(2 + n) 

+ 6(40 + 17n)) + (10a + 96) /r) + en((6 - 7an) cot 2 + 8a esc 2 0))) ] | ro , 

G(e\ ro ) = — l—jz . [ 2(10a + 76)/ 3 3 £ r 6 - ./ 2 fer 3 (20a(2 + n) + 66(6 + n) (A.9) 

2(18a + 176)/" ! er 4 

+ (26a + 7b)9nr cot 6) + 2n cot 9 (en(b + an) - r 2 + 2a/er 2 - e(-2 + n) (6 + an) esc 2 0) 
+ 2fr(-2en(2a + b - 2an) cot 9 + 9r(en(b - 7an) + r 2 - 4(2a + 6)/er 2 
+ en(-8a - 6 + 7an) esc 2 9)) } \ ro , 

H{9\ ro ,9\ r „) = ^- [ 8(6 - a(-2 + n))(-l + n)?icot 6» + (l/e)r((98a + 836)/ 3 6» 5 £ r 7 (A.10) 

128(a + 6)/r 4 

- 2(35a + 186)/ 2 6» 4 enr 5 cot 6> + 40r(en(-6 + 7an) - r 2 + 2er(/(20a(2 + n) 

+ 6(40 + 17n)) + (10a + 96) /V) + e(6 + a(8 - 7n))n esc 2 9) + A9{en 2 {-b + 7an) 

+ (-2 - n)r 2 + 2er(r(2/(3a(2 + n) + 6(5 + n)) + (2a + 6)/ (3) r) + /(2a(2 + n)(2 + 3n) 

+ 6(8 + n(15 + n)) + (28a + 196)07ir 2 cot 8)) + e(b + a(8 - 7n))?i 2 esc 2 9) 

+ 2f9 3 r 3 (en(b - 7an - 46n) - 3r 2 + er(/(-6o(34 + 7n) - 6(174 + 19n)) 

- 4(8a + 7b) fr) + en(b(-5 + in) + a(-8 + 7n)) esc 2 9) + 26> 2 r(-(210a + 197b)9f 2 er 4 
+ 2n cot 6»(en(-36 + an) + 2r 2 + er(/(66a + 376 + 4(6a + b)n) + (9a + 46) fr) 

+ e(3b - o(-4 + n))(-2 + n) esc 2 0))) ] | . 
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